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ABSTRACT 

We present a detailed comparison between the 2 — 10 keV hard X-ray and infrared (IR) lu- 
minosity function (LF) of active galactic nuclei (AGN). The composite X-ray to IR spectral 
energy distributions (SEDs) of AGN used for connecting the hard X-ray LF (HXLF) and IR 
LF (IRLF) are modeled with a simple but well tested torus model based on the radiative trans- 
fer and photoionization code CLOUDY. Four observational determinations of the evolution of 
2 — 10 keV HXLF and six evolution models of the obscured type-2 AGN fraction (fz) have 
been considered. The 8.0 and 15 [am LFs for the total, unobscured type-1 and obscured type-2 
AGN are predicted from the HXLFs, and then compared with the measurements currently 
available. We find that the IRLFs predicted from HXLFs tend to underestimate the number of 
the most IR-luminous AGN. This is independent of the choices of HXLF and fa, and even 
more obvious for the HXLFs recently measured. We show that the discrepancy between the 
HXLFs and IRLFs can be largely resolved when the anticorrelation between the UV to X-ray 
slope a ox and UV luminosity Ltjv is appropriately considered. We also discuss other possible 
explanations for the discrepancy, such as the missing population of Compton-thick AGN and 
possible contribution of star-formation in the host to the mid-IR. Meanwhile, we find that the 
HXLFs and IRLFs of AGN can be more consistent with each other if the obscuration mech- 
anisms of quasars and Seyferts are assumed to be different, corresponding to their different 
triggering and fueling mechanisms. More accurate measurements of the IRLFs of AGN, espe- 
cially that determined at smaller redshift bins and more accurately separated to that for type-1 
and type-2, are very helpful for clarifying these interesting issues. 

Key words: galaxies: active - galaxies: luminosity function - galaxies: formation - galaxies: 
evolution - infrared: galaxies - X-rays: galaxies 



1 INTRODUCTION 

Active galactic nucleus (AGN), compact regions at the center of 
active galaxies, releasing a great deal of energies in the form 
of radiation over the electromagnetic spectrum from radio, in- 
frared, optical, ultraviolet, X-ray to 7-ray, are now believed to 
be powered by accretion of mass into the super-massive black 
holes (SMBHs). In the local Universe, SMBHs are found to ex- 
ist at the center of most massive galaxies. There are good cor- 
relations betw een the mass of SMBHs and the proper t ies of 
host galaxies ( Hopkins et al. 2007; Kor mendv & Bender] 120091 ; 
Gultekin et al. 2009; Merloni et al. 2010) such as the veloc- 
ity dispersion JFerrarese & Merritfl [20 00; Gebhardt eta?] I2OO0I: 
Tremaine et alj|2002). mass JMagorrian et a l. 1998; H aring & Rixl 



20041 : lGrahaml2004l) or luminosity {Kprmendy & Richstonel 19951 : 
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iMarconi & Huntll2003l) of host bulge. On the other hand, the AGN 
activity and star formation are found to peak at a si milar red- 
shift and decline towards low redshift simulta neously dHopkinsI 
I2004J ; ISilverman et"aI1 120081 ; lAird et alj l20ld) . Meanwhile, the 
mass density of local SMBHs in galaxy center is found to be 
consistent wit h that accreted by AGN throughout the history of 
the Universe JMarconi etail |2004| ; iMerlonil I2004T) . These cor- 
relations strongly support the idea that the growth of SMBHs 
sho uld be coupl ed with the formation and evolution of galax- 
ies JCroton et al. 2006; Bower et al. 2006; Pi Matteo et al.ll2005L 
2008; Hop kins et al.1 120051. I2006T 20081) , although some authors 
(e.g. |Peng|2007l : |jahnke & MaccioluOllh argued a non-causal ori- 
gin of them. 

While the important role of SMBHs, and so AGN, playing 
in the formation and evolution of galaxies has been well estab- 
lished, detailed mechanisms about this process are still largely 
unknown. The luminosity functions (LF) of AGN, which de- 
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scribe the spacial density of AGN as a function of luminosity 
and redshift, is an important observable quantity for understand- 
ing the distribution and evolution of AGN. It constrains the ac- 
cretion history of SMBHs, and reveals the triggering and fuel- 
ing mechanism of AGN and their co-evolution with host galax- 
ies. An observational determination of the bolometric LFs of AGN 
require multi-wavelength observations spanning the whole wave- 
length range of electromagnetic spectrum and sampling large co- 
moving volume and luminosity range. So, in practice, the LFs of 
AGN are measured independently from different wavelength bands 
such as radio (e.g. Nagar et al. 2005), infrared (e.g. Babbedge et al. 
20061 : Brown et al. 2006; M atute et all2006t) . optical (e.g.lFan et all 



2001 



Wolfetal. 2003: 



Croom et al. 



20041; iRichards et al] 1200(1 



Bongiornoetalj 12007; Fontan oTeTalJ 120071: Ishankar & M athui 
2007b. soft X^ray (e.g. iMivaii et al] 12000, 2001; Sil verman et al 



2005al: lHasinger et alJ 120051). hard X-ray (e.g. lUedaet al]|20o1 
La Franca et alj|2005l:lsilverman et al. 2005bl 120081 ; I Yencho et al] 
20091; lEbrero et alj|2009l ; lAird et al]|2010h . or emission lines (e.g. 



Hao ct 



alj|2005l) . 



However, due to the different selection effect suf- 
fered by different bands, the LFs of AGN measured from different 
bands are not necessarily consistent with each other. 

Among various bands, X-ray, especially the hard X-ray band, 
is the most efficient for selecting AGN. Recently, the evolution of 
the hard X-ray LF (HXLF) of AGN from z ~ to 5 is found 
to be best described by a luminosity dependent density evolution 
(LDDE) model. According to this model, the spatial density of 
AGN with lower luminosity peaked at lower redshift than those 
with high luminosity, and t he faint-end slope of the LFs is flat- 
tened as red shift increased dUeda et alj|2003l : iBarger et al] 120051 : 
lHasinger et alj|2005t) . This kind of so-called 'cosmic downsizing' 
evolution trend of the AGN population has been further confirmed 
in rad io and optical bands dCirasuolo et alj|2005t iBongiorno et al] 
120071) . These results revealed a dramatically different evolutionary 
model for Seyfert galaxies and quasars, and imply very different 
triggering, fueling and accretion mechanisms for the two classes of 
AGN. 

Meanwhile, AGN are classified to two major classes accord- 
ing to their optical spectra. Type-1 AGN exhibits both broad per- 
mitted lines and narrow forbidden lines in their spectra, while type- 
2 AGN presents only th e narro w lines JKhachikian & Weedmar] 
1 1974ft . iRowan-Robinsonl d 19771) firstly put forward the idea that 
AGN are surrounded by dusty medium which absorbs their visi- 
ble and ultraviolet light and then re-emits them in the mid-IR. The 
extinction due to these obscuring medium is responsible for the 
distinction between type-1 and type-2 AGN. Latterly, this idea was 
develo ped into the so-ca ll ed unified model of AGN (Pier & Krolik 
2l; lAntonuccl Il993l: iMaiolino & Riekel Il995l : Would 1 19991 
Zhang & Wa ng 2006; Wang & Zhang 2007). In the model, the dif- 



ferences between different types of AGN can be explained by the 
anisotropically distributed obscuring mediums (often visualized as 
a geometrically and optically thick torus comprised of dust and 
molecular gas) surrounding a basic black hole-accretion disk sys- 
tem, while different lines-of-sight into and through these obscuring 
mediums result in the diverse observational properties of AGN pop- 
ulation. 

However, the obscuration of AGN by anisotropically dis- 
tributed gas/dust medium imply great systematic selection bias for 
understanding the properties and evolution of AGN. Moreover, 
obscuring mediums around AGN are recently found to be dis- 
tribute d in a much more complex manner than a simple compact 
torus (Risaliti et al. 2002; Kuraszkiewicz et al. 2003; Risaliti et al. 



minositv dUedaet al] 12 003; St effen et al] 120031 ; iHasingeil |2004|; 
Bargeretal. 2005; Simpson 2005) and redshift ( La Franca et al. 



2005; Ballanty neet al] l2006al : iTreister & Urrvl 120061 : lHasinge 



2009) 



1 
;er 



120051 ; iGoulding & Alexander 120091) . and may evolve with lu- 



2008; Ebre ro et al] 120091) . These results imply that the observa- 
tional properties of AGN may vary significantly from object to 
object. This complicates the understanding of the intrinsic char- 
acteristics of AGN and their correlations with the host galaxy. 
On the other hand, current state-of-a rt synthesis mod els of cos- 
mic X-ray background (CXRB) (e.g. iGilli et al]|2007]) show that 
a large population of heavily obscured Compton-thick AGN (with 
Nu ^ 10 24 cm -2 ) are required to fit the CXRB spectrum. This 
population of Compton-thick AGN can be missed by even the 
deep hard X-ray surveys since they are deeply buried by obscur- 
ing medium. 

Furthermore, recent results of iHasingen ( 1200 81) and 
ITreister et al] d20101) show that the fraction of absorbed AGN 
increases significantly with redshift to z ~ 2 — 3, accompanied 
with the cosmic co-evolution of star-formation and AGN activity. 
These results support the idea that the obscuration of AGN cannot 
simply come from an unevolving torus employed by traditional 
unified model of AGN. The obscuration mechanism of AGN 
with different triggering, fueling and accretion mechanisms may 
be different and associated with their co-evolution with ga laxies 
JDavies et aT1l2006l : lB"aiiantvne et alj2006al : lBallantvnel2008l) . 

So, detailed studies of the obscuring medium around AGN, 
such as their geometry, distri bution, composition, origin, and 
evolution, are ver y important (Zhang] [2004J; IWang et al] 120051 : 
iLiu & Zhandl201 ll) . The obscured or absorbed optical, ultraviolet, 
and X-ray radiation will be re-emitted in the infrared (IR). IR bands 
represent an important complement for understanding the proper- 
ties of obscuring medium around AGN and their co-evolution with 
host galaxies. With the existing IR space telescope such as Spitzer, 
Herschel and forthcoming James Webb Space Telescope (JWST), 
our observation and understanding of AGN from the IR band will 
be largely improved. Given the limitations suffered by X-ray ob- 
servations, it is important to study the LFs and obscuration of AGN 
together and test the conclusions about their evolution, which are 
mainly based on observations in the X-ray band, in the IR band. 

By using the spectral energy distributions (SEDs) modeled 
with a simple torus model, which is bas ed on the radiative trans- 
fer and photoionization code CLOUDY, iBallantvne et al] J2006bh 
can relate the X-ray and IR properties of AGN and explore the ef- 
fects of parameters about obscuring medium. They presented the 
mid-IR number counts and LFs for three evolution models of the 
$2 (equal to the covering factor under the unified model of AGN) 
that are constrained by the synthesis model of CXRB. The mid-IR 
number counts and LFs predicted from HXLF are in good agree- 
ment with direct IR observations, especially when assuming an in- 
ner radius (Ri n ) of 10 pc for the obscuring medium as expected if 
the obscuring material is connected to galactic-scale phenomenon. 
The mid-IR LFs of AGN are found to be a much better tool for 
determining the evolut ion of J 2 with z. 

Balla ntvne et al] d2006bl) presented the mid-IR LFs for to- 
tal AGN at d ifferent redshifts, b ut the observational mid-IR LFs 
of AGN (i.e. iBrown et alj|2006l) . which is used to be compared 
with, a re for type- 1 AGN only. After the work of IBallantvne et al] 
(2006b), some important improvements t o the measurement of 
HXLF of AGN have been presen ted (e.g. Isilverman et al . 2008; 
lEbrero et al . 2009; Aird et al. 2010). Furthermore, the actual evolu- 
tion model of AGN obscuration i s not necessa rily within the three 
models proposed by IBallantvne et al] ]2006al) , other possibilities 
need to be tested for more reasonable conclusions. 
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In this paper, we present a more detailed comparison between 
the HXLFs and mid-IR LFs of AGN, which are connected by the 
composite X-ray to IR SEDs modeled with a modified version o f 
the simple by well tested torus model of B allantvne et al.l 12006b). 
More observational determinations of the 2 — 10 keV HXLF and 
the evolution models of f% have been considered. The 8.0 and 
15 u.m LF for the total, unobscured type-1 and obscured type-2 
AGN are predicted from different combinations of HXLF and fa , 
and then compared with current I R ob servational results. Besides 
the measurement o f iBrown et al.l 120060, the 15 um L F given by 
iMatute et al] J2006I) and recent results of lFu et al.l J20ld) have been 
added for comparison. 

We begin in Section [2] by reviewing current understanding 
of AGN evolution from X-ray band. This include current obser- 
vational determination of the evolution of the HXLF of AGN in 
Section |2~T1 and the evolution of AGN obscuration in Section |2T2l 
The detailed procedures of modeling the composite X-ray to IR 
S EP of AGN, and our mo difications to the original torus model 
of iBallantvne et alj 12006b 5 ), are presented in Section [3] Section [4] 
presents the method used to compute the IRLFs of type- 1 , type-2 
and total AGN from different combinations of HXLF and fa. In 
Section|5] we present our results and compare them with measure- 
ments from direct mid-IR observations to seek conclusions about 
the evolution of LFs and obscuration of AGN from combined views 
of hard X-ray and mid-IR. We find that the mid-IR LFs predicted 
from HXLFs tend to underestimate the number of the most IR- 
luminous AGN, which is independent of the choices of HXLF and 
fa, and even more obvious for the HXLFs recently determined. In 
Section[6] we discuss explanations for this. Finally, a summary of 
this paper is presented in Section|7] 

Throughout this paper, we adopt a Hp = 70 km s _1 Mpc -1 , 
Qa = 0.7, and Q m = 0.3 (Spergel et al. 2003) cosmology. Minor 
differences in the cosmology have negligible effects on our conclu- 
sions. 



2 THE EVOLUTION OF AGN REVEALED FROM X-RAY 
BANDS 

2.1 The evolution of the HXLF of AGN 

Strong X-ray emission is a unique indication of an AGN activity 
at the center of galaxies. Deep X-ray surveys by Chandra and 
XMM-Newton, which have already resolved most of the 2 — 10 keV 
cosmic X-ray background (CXRB) into individual sources, found 
that most sources of CXRB are AGN. X-ray, especially hard X-ray 
with energy ^2 keV, are highly efficient for selecting AGN. Both 
the moderately obscured (A?h ^5 10 23 cm -2 ) and low-luminosity 
sources commonly missed by optical observations can be selected 
from hard X-ray. So, much more trustable evolution trends of the 
AGN can be revealed from hard X-ray observations. 

As mentioned above, HXLF of AGN is found to be best de- 
scribed by a LDDE model. However, the exact form of the evolu- 
tion is still under debate, especially at high redshifts. In this paper, 
we adopt the LDDE model given bv lUeda et alj 120031) . where the 
present-day HXLF is described as a smoothly-connected double 
power-law form: 
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Figure 1. HXLF of AGN a t z = 0.5, 1.5,2.0 and 2.5 as given by the 
LDDE model of llJeda et alj 120031) (red solid line, noted as U 03LDDE) 
lEbrero et alj J2009l) (green dashed line, noted as E09LDDE), lAird et alj 
bOld) (blue short d ashed line, noted as A10LDDE) and the LADE model 
of lAird et al] 12010) (purple dot line, noted as A10LADE), respectively. 



evolution of LFs is given by 

d$(Lx,z) _ d$(L x ,0) 



e(z,L> 



(2) 



dlogLx dlogLx 

where the evolution term is given by 

( (l + z) pl [z<z c (L x )] 

e(z,L x ) = \ (3) 

I e (^)[irifer 2 [z>z B (L x )]. 

The cutoff redshift z c , with a dependence on the luminosity starting 
from a characteristic luminosity L a , are given by a power law of 
L x : 



Z c {Ly 



(L X > L a ) 



z'A^) (L x <L a ) 



(4) 



where 7I is faint-end slope, 72 is the bright-end slope, L, is the 
characteristic break luminosity and A is a normalisation factor. The 



where a measures the strength of the dependence of z c with lumi- 
nosity. 

Recently. lEbrero et alj J2009I) re-measured the HXLF of AGN 
by us ing the XMM-Newton Medium Survey (XMS, Barcons et alj 
2007) and other highly complete deeper and shallower surveys to 
assemble an overall sample of ~ 450 identified AGN in the 2 — 10 
keV band, which is one of the largest and most complete sample up 
to date. lAird et alj 120101) presented a new observational determi- 
nation of the evolution of the 2 — 10 keV HXLF of AGN by using 
data from many surveys including 2 Ms Chandra Deep Fields and 
the AEGIS-X 200 ks survey. These, combined with a sophisticated 
Bayesian methodology, allow them to do a more accurate measure- 
ment of the evolution of the faint end of the HXLF. They found that 
the evolution of HXLF are best described by a so called luminosity 
and density evolution (LADE) model, rather than the LDDE model. 
The LADE model is a modified Pure luminosity evolution (PLE) 
model. According to the PLE model (Ueda et al. 2003), the evolu- 
tion of HXLF with redshift are described by allowing the charac- 
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teristic break luminosity L« in the present-day HXLF as given by 
Eq. (0 to evolve as 



log L,(z) = logLo -log 



l + z c 
1 + z 



l + z c 
1 + z 



(5) 



where the parameter z c controls the transition from the strong low- 
z evolution to the high-z form. The LADE model are constructed 
by additionally allowing for overall decreasing density evolution 
with redshift, i.e. allowing the normalization constant A in the 
present-day HXLF as given by Eq. (T) to evolve as 



log A(z] =logA +d(l + z) 



(6) 



where d is an additional parameter describing the overall density 
decreasing. 

In Fig.Q] we show the HXLF of AGN at z = 0.5 1.5, 2.0 and 
2.5 as g iven by the L DDE model o flUeda et all fc003l) . lEbrero et al 



2009), lAirdetalJ ( feoioh . and the LADE model of lAird et al 



20 1 0) . The LADE modeling of the HXLF retains the same shape at 



all redshifts, but undergoes strong luminosity evolution out to z ~ 
1. Meanwhile, the HXLF undergoes overall negative density evo- 
lu tion with in creasing redshift. Fig.[T]clearly shows that the HXLF 
of lAird et al.l 120101) is different from those of the others at high 
luminosity. 

Different HXLFs have very different implications for the AGN 
populations, such as their lifetimes, duty cycles, fueling, triggering 
and evolution. Further complemental views from other wavelength 
bands, such as IR, are important for a fully understanding of the 
evolution of AGN populations. In this paper, the four modelings of 
HLXF mentioned above have been used to predict the evolution of 
mid-IR LFs of type- 1 and/or type-2 AGN, respectively. 



2.2 The evolution of the obscuration of AGN 



Since there is a good correspondence between the AGN with 
Nh z2 10 22 cm" 2 and those optically identified as being of type-2 
JTozzi et al. 2006), type-2 AGN is commonly defined as those with 
absorbing column densities Nn ^ 10 22 cm -2 in the X-ray band. 
According to the unified model of AGN, fy approximately equal to 
the covering factor of the gas with Nu ^ 10 22 cm -2 around the 
AGN. 

By using a pop ulation synthesis model of CXRB, 
iBallantvne et alj <2006al) constrained the evolution of f 2 as a 
function of both z and Lx- They presented three parameterizations 
for the evolution of f 2 (logLx , z) that could fit the observed shape 
of the CXRB and X-ray number counts of AGN. The first one 
(shown in Fig. [2] and noted as 'f2_l'), with a moderate redshift 
evolution, is given as: 



/ 2 = Ki(l + z)°- 3 (log-L s 



(V) 



where Ki is a constant defined by /2(logLx = 41.5, z — 0) = 
0.8, which is based on observations in the local Universe. The sec- 
ond one (shown in Fig. [2] and noted as 'f2_2') with a more rapid 
redshift evolution, is given as: 



/2 = K 2 (l + 2f 9 (lo gj Lxr 



(8) 



where K2 is based on the Sloan Digital Sky Sur vey (SPSS) mea - 
surement of / 2 (logLx = 41.5, z= 0) = 0.5 bv lHao et~aIU2005l) . 
In the above two cases, the z evolution is halted at z — 1, be- 
cause there is no constraint on f 2 at higher redshifts. The last one 
(shown in Fig. [2] and noted as 'f2_3'), which is considered as a 



null-hypothesis, assumes that f% does not evolve with redshift, and 
given as: 



/ 2 = K 3 cos 



logLx - 41.5 
9/7 



(9) 



where K3 is determined by /2 (logLx = 41.5, z — 0) = 0.8. 

The fraction of type-2 AGN can also be measured directly 
from observations in different bands. For example, in the op- 
tical, AGN can be selected using their high ionizat ion lines to 



construct the stand a rd diagnostic dia grams (Baldwin et al. 1981; 



lKewlevetal1l200ll: llCauffman n et alj[2003t [Kewley et al 



Furthermore, the ratio of narrow-line and broad-line AGN can be 
measured as a function of 2 and the luminosity of emission lines 
(such as [O III] 5007A line), which can be used as AGN power 
indicators. However, significant limitations of optical s election and 
classification of AGN have been noti ced (e.g. iMoran et alj|2002t 
iNetzer et alJuOOfl ; iRigbv et alj |2006). The nuclear emission can 
be obscured by the torus and/or outshone by the host-galaxy light. 
Alternatively, AGN can be efficiently selected in the X-ray, and 
the X-ray luminous AGN can be classified to absorbed and unab- 
sorbed according to their absorbing column densities logA^n < 22 
or > 22. X-ray selection of AGN is suffered by the limited sen- 
sitivity of telescopes, which are only sensitive at ^ 10 keV. So, a 
significant fraction of absorbed objects, especially the large number 
of Compton-thick AGN with logA^H > 24 predicted by the popu- 
lation synthesis model of CXRB, may be missed by current hard 
X-ray selection of AGN. 

The combination of X-ray and optical criteria is a much more 
robust method for the selection and classification of AGN. Re- 
cently, |HasingerJ {2008) have presented a new determination of the 
fraction of absorbed sources as a function of X-ray luminosity and 
redshift from a sample of 1290 AGN. They are selected in the 2 — 10 
keV band from different flux-limited surveys with very high opti- 
cal identification completeness, and grouped into type-1 and type- 
2 according to their optical spectroscopic classification and X-ray 
absorption properties. So, the evolution of AGN absorption with lu- 
minosity and redshift is determined with higher statistical accuracy 
and smaller systematic errors than previous results. The absorbed 
fraction is found to decrease strongly with X-ray luminosity, and 
can be represented by an almost linear decrease with a slope of 
0.281 ± 0.016. Meanwhile, it increase significantly with redshift 
as ~ (1 + 2 ) 062±011 from z = to z ~ 2. On the other hand, 
the evolution of the absorbed AGN fraction over the whole redshift 
from z = to z ~ 5 can also be described as ~ (1 + z )°- 48 ± 008 5 
or ~ (1 + 2;) - 38 ± 009 w hen data with crude redshifts are excluded. 

These findings may have important consequences for the 
broader context of AGN and galaxy co-evolution. According to the 
results of Hasinger (2008), we have constructed three new evolu- 
tion models of AGN obscuration, which are expressed as 

sO.62 



h 



-0.281(logL x - 43.75) + 0.279(1 + z) 



f 2 = -0.281(logL x - 43.75) + 0.308(1 + z) 



h 



-0.281(logL x - 43.75) + 0.309(1 + z) 



(10) 



(11) 



(12) 



and noted as 'f2_4', 'f2_5', and 'f2_6\ respectively (shown in Fig. 
O. Due to the simple linear dependence on luminosity, the type-2 
AGN fraction will quickly bec ome zero as lumin osity increasing. 
According to recent results of iBrusaet all J20101) . the fraction of 
the obscured AGN population at the highest (Lx > 10 44 erg s _1 ) 
X-ray luminosity is ~ 15% — 30%. So, we have set a lower limit 
of 0.15 for the evolution of f 2 (logLx, z) to stand for a flattening 
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Figure 2. The six evolution models of /2(log£xiZ) used in this paper are shown at logLx=41.5 (red solid line), 43.0 (green dashed line), 44.5 (blue 
short dashed line) a nd 46.0 (purple dot line), respectively. The first three, which are constrained by CXRB spectrum and X-ray number count as given 
by Ballantv ne et alj <2006al) . are noted as '£2_1', 'f2_2', and 'f2_3', respectively. The last three, which are constructed according to recent measurement of 
lHasingen 120081) . are noted as 'f2_4', 'f2_5\ and 'f2_6', respectively.(see text for more detailed explanations for these evolution models.) 



of the decline at the highest luminosities as expected. Naturally, an 
upper limit of 1 is forced at all cases. 

Except for /2(logLx, z), additional assumptions are needed 
to determine a specific distribution of Ah- The exact distribu- 
tion of TVh is un known except for local bright Seyfert 2s (e.g. 
iRisaliti et al . 1999), but the covering factor is a useful parameter 
for its theoretical description. Here, we use a simple as s umptio n 
about the distribution of Ah following iBallantvne et alj {2006b). 
In the 'simple Nh distribution', ten values of Ah are considered: 



log (Ah /cm" 



20, 20.5, . . . , 24.0, 24.5, and a type-1 AGN 



is assumed to have an equal probability p of being absorbed by 
columns with logA/H < 22. Likewise, a type-2 AGN has an equal 
chance of being absorbed by columns with logA^H ^ 22: 



log(Ar H /cm 



20.0, 



22.0, 



,21.5, 



,24.5, 



l-/2Qog-f-X, z ) 



J 2 (logLx,z) 
6.0 



(13) 

Since /2(logLx,z) depends on logLx and z, the distribution of 
obscuring medium around AGN evolves with both logLx and z. It 
is worth noting that this simple assumption about the distribution 
obscuring medium with different Ah is only used to construct Ah- 
averaged SEDs (described in the next Section). We do not expect 
it to give a correct fraction of Compton-thick AGN. In fact, there 
are AGN with e stimated column logAfa > 25 (e.g. NGC 1068; 
uMatt et al .1119971) . The inclusion of very Compton-thick obscur- 
ing medium dramatically increase the computation time of SEDs. 
However, they have important effects mainly in the far-IR, but only 
ignorable effects in the mid-IR which we are mostly interested in 
currently. 



3 MODELLING THE SPECTRAL ENERGY 
DISTRIBUTION OF AGN 

To predict the IR properties of AGN, we must know the relation be- 
tween IR and X-ray luminosity of AGN. This can be given by the 
SEDs of AGN with different X-ray luminosities. The SEDs of AGN 
can be obtained from observations or theoretical calculations. The 
observational SEDs have the advantage of being based on obser- 
vations of real AGN. However, the number of observed objects is 
limited, and they only cover a narrow range of luminosity, redshift 
and wavelength. Alternatively, we can use theoretical dusty torus 
emission models, which include a detailed radiative transfer calcu- 
lation, to compute the expected IR SED for a given X-ray luminos- 
ity. However, most radiative transfer calculation for IR dust emis- 
sion of AGN do not incl ude detailed considerations of gas and its 
interaction with dust (e.g. lTreister et al . 2004, 2006; N enkova et alj 



2008a. b). The gas and dust are expected to be interacting with each 
other, and gas is responsible for the absorption of X-ray. So, for a 
reasonable connection of the X-ray and IR properties of AGN, gas 
and their interaction with dust must be considered. 

Here, the calculation of AGN SEDs is performed by using the 
photoionization code CLOUDY v. 07.0 2.01 JFerland e t al. 1998), 
following a procedure similar to that of IBallantvne et alj |2006b), 
but with some simplifications. In CLOUDY, the atomic gas physics 
along with the detailed dust radiation physics, such as poly- 
cyclic aromatic hydrocarbon (PAH) emission and emission from 
very small grains, have been self-consistently considered. In addi- 
tion, many important physical properties of the obscuring medium 
around AGN, such as its distance from the central engine, gas den- 
sity, distribution and gas/dust ratio, can be varied freely, and so ex- 
plored extensively. CLOUDY is a one-dimensional radiative trans- 
fer code, and the methods we employed to model the SED s of AGN 
are less sophisticated than those used by [Treister et al.| d2006l) or 
iNenkova etal] fc008al) . However, IBallantvne et alj J2006bl) showed 
that the SEDs, when averaged over a A^h distribution, have very 
similar properties to the ensemble of AGN found in the deep sur- 
veys of Chandra, XMM-Newton and Spitzer. 



3.1 Construction of CLOUDY model 

To construct a CLOUDY model, three ingredients must be speci- 
fied. Firstly, the shape and intensity of the radiation source, which 
define the incident continuum, must be set. The intrinsic spectrum 
of AGN is described by a multi-component continuum typical for 
AGN, which extend from 100 keV to > 1000 ^.m. Specifically, 
the 'Big Bump' component, peaking at 1 Ryd, is a rising power 
law with a high-energy exponential cutoff and parameterized by 
the temperature of the bump. The big blue bump temperature is set 
to be a typical value of 10°K. The X-ray to UV ratio a ox , which is 
defined by 



log(L2kcv/£ 2 500A) 
log02keV/V 2 500A) 



0.3838 log 



J 2500A 



(14) 



have an important effect on the result ing X-ray to IR ratio. Espe- 
cially , there are evidences (e.g. iSteffen et alj|200q ; iHopkins et alj 
120071 ; IVagnetti et alj|20101) that this parameter may be anticorre- 
lated with the UV luminosity of AGN. To explore the effects of 
this important parameter, we set q ox to be a constant value of — 1.5, 
— 1.4 and — 1.3, respectively. We a lso tested the q ox — Luv relation 
presented in lHopkins et alj 020071) . which is given by 



aox = -0.107 log( 



ergs -1 Hz" 



+ 1.739, 



(15) 
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and determined specifically for unobscured (type-1) quasars. This 
results in a luminosity-dependent shape of the input SED from 
AGN center. The low-energy slope of the Big Bump continuum 
a uv is set to be the default value of —0.5. The X-ray photon index 
is assumed to be r = 1.9, and so the energy index a x = 1 — V = 
—0.9. The full continuum is the sum of two components as given 
by 



u 



exp(-hv/KT B B) exp(-KTm/hv) + av° 



(16) 



where Tbb is the temperature of Big Bump and the coefficient a 
is adjusted to produce the correct <a ox for the case where the Big 
Bump does not contribute to the emission at 2 keV. The Big Bump 
component is assumed to have an IR exponential cutoff at KTir. = 
O.OlRyd (lRyd ~ 13.6 eV)|j Finally, this spectrum is scaled to 
have a luminosity of L 2500 ^(erg s _1 Hz -1 ). 

The second ingredient of a CLOUDY model is the chemical 
composition of the obscuring medium. A gaseous element abun- 
dance similar to that of Orion Nebula is assumed. The size distribu- 
tions and abundances of graphitic, silicate and PAHs grains are also 
set to be similar to that of Orion Nebula. The obscuring medium is 
assumed to distribute uniformly and have a constant hydrogen den- 
sity nH of 10 4 cm -3 . 

The last ingredient of a CLOUDY model is the geometry of 
the obscuring medium. Here, the obscuring medium is assumed 
to be Ri n pc away from the center and with a column density of 
Ah. On the other hand , to be consistent with the unified model, 
iBallantvne et alj ( 120060) have set the covering factor of the obscur- 
ing medium to fa when Ah ^ 10 22 cm -2 or 1 — fa otherwise, in 
CLOUDY model. Since fa depends on both luminosity and red- 
shift, the CLOUDY simulation needs to be done for each luminos- 
ity and redshift, respectively. This would result in a great number of 
CLOUDY models. However, in CLOUDY models the covering fac- 
tor only has second-order effects on the spectrum through changes 
in the transport of the diffuse emission. So, we just use the de- 
fault geometric covering factor of unity (the shell fully covers the 
continuum source) but a radiative covering factor of zero, i.e. an 
open geometry is assumed, and the reflected radiation can be ob- 
tained as well. The effects of covering factor on the diffuse and 
reflected emissions are considered after the CLOUDY simulation 
as described int the next section. 



3.2 CLOUDY model grids and construction of AGN SEDs 

The CLOUDY models are built for i 2 50oA( er S s_1 Hz _1 )[ 2 Jfrom 
27 to 34 (in steps of 0.2 5), and loR(A H /cm 2 ) f rom 20.0 to 24 (in 
steps of 0.5). Following Bal lantvne et al] I l2006bl) . we firstly set R ln 
to be 10 pc. However, we found that the temperatures of grains will 
be much higher than their sublimation temperatures at the high- 
luminosity end if a constant Ri n of 10 pc is assumed. We have 
practically found a luminosity-dependent _Rj n , which is given by 



r t /£„(2500A) ll/2 

R ™ = 10 * I 7^46 1 P°' 



(17) 



to fix this problem. The CLOUDY models are also built by as- 
suming luminosity-dependent R[ n for comparison. Finally, as men- 
tioned above, the CLOUDY models are built for four choices of 
Qox, respectively. 

1 See CLOUDY document for more detailed explanations for the construc- 
tion of this AGN spectrum. 

To explore the effects of Ljjv-dependent R ln and o Q x on the resulting 



For each CLOUDY model, three kinds of SEDs are predicted: 
the attenuated incident continuum, diffuse continuum and reflected 
continuum. The SEDs with different Ah correspond to observa- 
tions from different direction. However, according to the unified 
model of AGN, obscuring medium with all values of column den- 
sity Ah simultaneously exist around AGN. On the other hand, 
there are evidences that obscured an d unobscured AGN present 
more similar Lmir/Lx ratios (e.g . Alonso-Herrero et ak[ 2001c 
iKrabbe et alJl200ll : iLutz et alJlJX)! iHorst et alj|200d l2008n than 
that predicted by traditional torus models assumin g a smooth dis- 
tribution of dusty obscuring medi um. Recent works {Nenkova et al.1 
l2008allbt lH6nig et alj|20"lfl iHonig & Kishimotoll2010l) show that 
the distribution of dusty obscuring medium is clumpy rather than 
contiguous. So, when we observe an AGN from one direction, both 
diffuse and reflected emission from all Ah can be observed, in 
addition to the attenuated incident emission through an obscuring 
medium with a column density Ah of this direction. For this rea- 
son, we have made a modificatio n to the original torus emission 
model of lBallantvne et al.1 J2006bl) . The SED of an AGN with col- 
umn density Ah is constructed by adding the diffuse and reflected 
emission averaged over all 10 models with different Ah to the at- 
tenuated incident emission through a particular Ah. The weights 
are given by the probability distribution of the column densities, 
which is a function of fa(logLx, Z) (or covering factor) as dis- 
cussed in Section 12.21 The SEDs constructed this way is called 
'unified SEDs'. 

Finally, the 'unified SEDs' undergo an average over Ah again 
to produce the 'Ah -averaged SEDs', which will be used to pre- 
dict the IRLFs of AGN later. Here, three types of Ah -averaged 
SEDs are constructed. The type-1 SED is an average of the 'uni- 
fied SEDs' with 10 200 cm~ 2 < N H < 10 22 cm~ 2 . The type- 
2 SED is an average of the 'unified SEDs' with 10 22 cm~ 2 ^ 
Nh ^ 10 24 ' 5 cm -2 . The average SED is an average of the 
'unified SEDs' with 10 200 cm- 2 sC N H < 10 245 cm -2 . As 
an example, Fig [3] shows the rest-frame SEDs taken from the 
'f2_l' evolutionary grid (eq. |7J for a Seyfert-like AGN (Lx = 
10 4354 erg s-\ 2 = 0.7, fa = 0.7484) and quasar-like AGN 
(Lx = 10 46M erg s" 1 , z = 1.4, fa = 0.5705), respectively. 



3.3 Testing model SEDs of AGN 

The metho d we have used to model the SEDs of AGN is similar to 
that of Ball antvne et al.1 J2006bh . They have been extensively tested 
this method against large samples of AGNs. Here, we present two 
additional tests that are more directly related to the goal of this 
paper, i.e. prediction of IRLFs from HXLF. For this goal, the most 
important thing i s correct X-ray to IR relation. 

Recently. iGandhi et al. (2009) found a strong mid-infrared:X- 
ray (logALA(12.3^m)-logZ/2_iukcv) luminosity correlation for a 
sample of local Seyferts, the cores of which have been resolved in 
the mid-IR. The relation is given by 



logAL A (12.3^m) 



-4.37 + 1.106 logLa-iotov, 



(18) 



and is found to be valid in a w ide range of luminosit y and may 
extend into the quasar regime. llVIullanev et al.1 ( 1201 ID converted 
this correlation to that between logLiR and logL2_iokcV, which 
is given by 



log- 



SEDs, L. 



2o()0A 



instead of Lx is used to define the luminosity of input SED. 



■ln43 =1 =0.53 + 1.11 log -—r — . (19) 

10 4,3 erg s L 10 4,i erg s i 

In Fig [4] the logLiR,-logI/2-iokeV relation computed from 
our model SEDs of AGN by using different choices of q ox 
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Figure 3. Top: Rest-frame SEDs for a Seyfert-like AGN (L x = 
ig43.5 er g s — 1) at 2 = 0.7 and obscured by dusty medium with an in- 
ner radius 10 pc, hydrogen density njj = 10 4 cm -3 and covering factor 
f<2 = 0.7484. The type-1 SED is shown in red, the type-2 SED is shown 
in green, while the average SED is shown in blue. Bottom: As top, but for 
a quasar-like AGN (Lx = 10 46,5 erg s — 1 ) at z = 1.4 and obscured by 
dusty medium with a covering factor f2 = 0.5705. These spectrums are 
taken from the 'f2_l' evolutionary grid (eq.|7J. 



and -Ri n are tested a gainst the observational relation given by 
iMullanev et alj d201 lh . As shown in the figure, th e result obtained 
for q ox = —1.4 and Ri n — 10 pc (as in iBallantvne et al.1 
(2006b)) significantly deviates from the nearly linear relation given 
by Mullanev et al. (201 1). This problem has also been noticed by 
iDraper & Ballantvnd ( 1201 II) . We found that a more linear logLiR- 
logZ/2-iokcV relation can be obtained if R ln decreases with Ltjv 
(as given by Eg. 117b. With the t ypical value of Qox = —1-4, this 
lead to a result similar to that of lMullanev et al.l J20111) . especially 
at the high-luminosity range. At the low-luminosity range, a larger 
Qox seems required. H owever, when the a x-£uv relation given by 
iHopkins et al.1 ( 120071) (as given by Eq. 115b is assumed, a too steep 
relation is obtained. So, these results support the anticorrelation be- 
tween q ox and Ltj v found by other independent observations (e.g. 
ISteffen et aT]|2006l : IVagnetti et"aT]|2010l) . but imply a more flat re- 
lation. 
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Figure 4. Test of the logLiR-logL2_ lrjkeV relation computed from our 
model SEDs of AGN against the observational relation (black-solid line) 
given by Mullanev et al. 12011). The results obtained by using different 
choices of a x and R ln are presented, respectively. The Ltjv -dependent 
-Ri n is given by Eq. 1171 while the Luv-dependent Qox is given by Eq. l 151 



4 CONNECTING X-RAY AND IR 

The IR is less affected by the selection effects due to obscuration 
suffered by optical and X-ray, while X-ray is currently the most ef- 
ficient for selecting AGN to high redshifts. Connecting X-ray and 
IR can provide a more clear view on the evolution of AGN popula- 
tions. If the evolution of AGN shown in HXLF and /2, which are 
revealed mainly from X-ray observations, are intrinsic, they should 
be shown somehow in the directly observed IRLFs of AGN. 



4.1 The mid-IR LFs of type-1 and/or type-2 AGN 

Since the HXLF tells us how the number density of AGN per incre- 
ment of logix changes with z and Lx, the following expression 
can be used to relate HXLF to IRLFs d$/d(log^L„): 

d$ = d$ d(logLx) 
d(logKLv) d(logLx) d(logz/L„) ' 

Here, d$/d(logL x ) is the HXLF of AGN given in Sectional] and 
L„ is the luminosity at a given wavelength. 

In Section [3] we have obtained the SEDs spanning from X- 
ray all the way to IR for AGN with different luminosities and red- 
shifts. So, the dependence of IR luminosity on hard X-ray luminos- 
ity, as described by d(logLx)/d(logz/L„), can be obtained from 
the SEDs easily. For predicting the IRLFs for total AGN (type-1 + 
type-2), we use the average SED presented in Section[3] The IRLF 
of AGN is not an integrated quantity, and so much more sensitive 
to the evolution trends in the HXLF and the /2 than the cumula- 
tive number count distribution and background spectra intensity of 
AGN. Using the IRLFs for total AGN has the advantage of being 
independent of the methods used to do a further classification of 
AGN, such as detailed optical emission line spectra, or an accurate 
measure ment of X-ray absorbing column densities TVh. For these 
reasons, Ballantv ne et al.1 ([2p06b) suggested to use IRLFs for total 
AGN to distinguish differe nt evolution models of A GN obscura- 
tion. However, the results of IBallantvne et al.1 J2006bl) showed that 
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the IRLFs is not very sensitive to the evolution model of AGN ob- 
scuration unless at much longer IR wavelengths where the contam- 
inant from star formation is important. 

The separated IRLFs for type-1 and type-2 AGN respectively 
are expected to be much more sensitive to the overall evolution of 
AGN spatial density and obscuration, although a detailed classifica- 
tion are required. The classification of AGN into type-1 and type-2 
involves with the problem of consistency between different classifi- 
cation methods. However, this may be a possible key to the problem 
of the evolution of AGN obscuration, since inconsistent classifica- 
tion methods will directly result in very different conclusions for 
the evolution of AGN obscuration. By separating the IRLFs to that 
for type-1 and type-2 AGN, the intrinsic evolution of AGN and the 
variations just resulted from the evolution of AGN obscuration, can 
be investigated in more detail and possibly clarified. Furthermore, 
by considering the IRLFs for type- 1 and type-2 AGN separately, the 
modeling of AGN SEDs can be further constrained. So, it would be 
more fruitful to separate the IRLFs of AGN into that for type-1 and 
type-2, respectively. This may provide a more useful tool to explore 
the properties of obscuring medium around different types of AGN. 

The separated IRLFs of type-1 and type-2 AGN are given as, 



d$i 



(1* 



d(\ogvL v ) d(logLx) 
and 

d$ 2 d$ 



(l-/a(logLx,«)) 



d(logLx) 
d(\oguL v ) 



d(logi/L„) d(logL> 



-/ 2 (logL x ,z) 



d(logLx) 
d(loguL„) 



(21) 



(22) 



where /a (logLx, z) is the fraction of type-2 AGN. We use the type- 
1 SED presented in Section[3]to predict the IRLFs of type-1 AGN, 
while using the type-2 SED to predict the IRLFs of type-2 AGN. 



5 RESULTS 

In this section, we present the predicted mid-IR LFs for total, type- 
1 and type-2 AGN, respectively. The SEDs of AGN used to obtain 
the X-ray to IR luminosity relation are computed by assuming a 
constant a ox = —1.4, and Luv -dependent R[ n as described by 
Eq.[T7] We leave the discussion of Luv-dependent a ox in Section 
[6] The results for all combinations of HXLF as given in Section l2~T1 
and different evolution models of AGN obscuration as given in Sec- 
tion |2J2] are presented and then compared with the measurements 
of mid-IR LFs of AGN currently available, respectively. Since we 
are mainly interested in finding out much obvious trends, a simple 
qualitative comparison by eye rather than a much detailed fittingQ 
is taken here. 



5. 1 The mid-IR LFs of total AGN 

In Figs.|5]and|6] the predicted rest-frame 8.0 and 15 tim LF for to- 
tal AGN are shown. In each panel, the results are predicted from an 
evolution model of HXLF and six evolution models of AGN obscu- 
ration as discussed in Section |2J2] The observational results used 
for comparison are from iFu et alj i201fj) . They used high-quality 
Spitzer 7 — 38 [im spectra to cleanly separate star formation and 
AGN in individual galaxies for a 24 (im flux-limited sample of 
galaxies at 2 ~ 0.7, and decomposed the mid-IR LFs between star 
formation and AGN. 

J Much more careful considerations of the covariance between points or 
systemic errors are not included as well. 



As can be clearly seen in Figs. [5] and [6] our results agree 
with that of Fu et al. (2010), Matute et al. (2006) and Hopkins et al. 
(2007) reasonably. These general agreements show that the meth- 
ods we have used to model the SEDs of AGN and to predict corre- 
sponding mid-IR LFs from HXLF are basically reasonable. Specif- 
ically, different evolution models of HXLF give very similar results 
at 8.0 /im and 15 (im. However, different evolution models of AGN 
obscuration are distinguishable at 15 /im, while not at 8.0 /im. As 
shown in Fig [6] the results at 15 /x m are divided into two groups, 
corresponding to using models from lBallantvne et alj (2006a) and 
using models constructed according to recent results of iHasingea 
(2008), respectively. The results predicted by using the e volution 



models of AGN obscuration from Ballantyne et al. 



pr, 
B 



better agreement with the measurements of lFu et alj J20101) and the 



2006a) are in a 



results from other authors, especially at the relative higher lumi- 
nosities. It seems that the e volution of AGN obscur ation are better 
described by the models of lBallantvne et alj d2006al) at th e redshift 
and luminosity ranges covered by the measurements of IFu et alj 
J2010h . i.e. 2 < 1 and uL„(8.0 /jm, 15 /jm) < 1O 12 L . 



5.2 The mid-IR LFs of type-1 AGN and type-2 AGN 

As mentioned above, it is more fruitful to separate the mid-IR LFs 
to that for type- 1 and type-2 AGN, respectively. Here, we present 
the mid-IR LFs for type-1 and type-2 AGN and then compare them 
with the mid-IR observational results of iBrown et al. (2006) and 
iMatute et all J2006t) . respectively. 



5.2. 1 The 8.0 \im LF of type-1 AGN 

From a sample consisting of 292 24 p_m sources brighter than 1 mjy 
selected from Spitzer MIPS survey, Brown et al. (2006) have deter- 
mined the rest-frame 8.0 \im LF for type-1 quasars w i th 1 < z < 5 
and 1.5 < z < 2.5, respectively. iBallantvne et alj (12006b) used 
these results (in their Fig. 13), but compared them with the pre- 
dicted mid-IR LFs for total AGN. Despite this, they found that the 
predicted and measured LFs s how a surprising agreement. As sug- 
gested bv lBrownet alj ( 120061) ■ if the fraction of obscured quasars 
decreases rapidly with increasing luminosity, the type- 1 quasar LF 
given by them would appropriate the LF of all quasars at the high- 
est luminosities. However, if there are indeed very few type-2 AGN 
at very high luminosities, and the type-1 AGN LF provide good ap- 
proximation of the total LF at the high luminosities, then this would 
be an important constraint for the evolution of AGN obscuration 
at high luminosities. T his important info r mation was not fully uti- 
lized in the method of lBallantvne et alj |2006b). So, it would be 
more reasonable and fruitful to predict the mid-IR LFs for only 
type- 1 AGN from different evolution models of HXLF and obscu- 
ration of AGN, an d then compare them with the measurements of 
lBrownetai]fc006h . 

In Fig. [7] we present the predicted rest-frame 8.0 (im LF of 
type-1 AGN at z = 1.5, 2.0, and 2.5, and compare them with 
the measurements of IBrown et alj d2006l) . As expected, the results 
predicted from different evolution models of AGN obscuration are 
more distinguishable when the rest-frame 8.0 ^m LF for only type- 
1 AGN, instead of total AGN, are used. The rest-frame 8.0 p.m LF 
of type- 1 AGN predicted from different choices of H XLF tend to 
undere stimate the number of AGN as measured by IBrown et alj 
(2006), and no matter which evolution m odel of AGN obscuration 
is used. Surprisingly, using the HXLF of lUeda et alj J2003|) resul t 
in a better agreement with the measurements o flBrown et alj (2006) 
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Figure 5. Rest-f rame 8.0 um LF for total AGN at z = 0.7 as predict ed from the LDDE modeling of HXLF of lUeda et all <2003l) , lEbrero et al.l J2009I) . 
lAirdet alj<2010h . and the LADE modeling of HXLF of Aird et al. 12010), respectively. At each panel, the results for six evolution models of /i are presented. 
The blue lines show the results for the three evolution models of J2 given by Ballantyne et al. 1 2006a), who constrained the evolution of f2 by fitting to 
the shape of CXRB spectrum and X-ray number counts. The short dashed-doted, long dashed-doted, and doted lines sh ow the results fo r the 'f2_l\ 'f2_2', 
and 'f2_3' evolution models, respectively. The red lines show the results for the three evolution models of J2 given by Hasinaer 12008) from direct X-ray 
observations of the evolution of type-2 AGN fraction. Here, the short-dashed, long-dashed, a nd solid l i nes sh ow the results for the 'f2_4', 'f2_5\ and 'f2_6' 
evolution models, respectively. The data points ar e the 8.0 um IRS-dec omposed A GN LF of iFu et al] 1201011 at z ~ 0.7. The green dot-dashed lines show 
the obscuration-cor rected AGN bolometric LF of iHopkins et al.l J2007I) (taken from | Fu et al J 1201(1) ). The purple solid lines show the LF for total AGN of 
iMatute et all 120061) as combined and converted to 8.0 urn at z ~ 0.7 bv lFu et alj fcOld) . 
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Figure 6. Similar to Fig. [3 but for the 15 u.m LF. The AGN obscuration evolution models of Ballantyne et al. 12006b) give better agreements with the 
measurements of Fu et al. (2010), Matute et al. 12006), and the results of Hopkins et al. 12007). 



than using all the more recent HXLF measurements. The results 
shown here may have, to some extent, confirmed the credibility of 
widely used results of Ueda et al. (2003). However, this may just 
be a coincidence, since more recent observational determination of 
HXLF are generally expected to be more accurate. More reasonable 
conclusions can only be obtained from comparisons with other in- 
dependent measurements. 

5.2.2 The 15 \un LF oftype-l/type-2 AGN 

From a samp l e of A GN selected at 15 (i.m (ISO) and 12 \xm (IRAS), 
IMatute et alj d200q) measured the rest-frame 15 (i.m LF of type-1 
and type-2 AGN, which are classified based on their optical spectra, 
separately. In Figs. [8] and [9] we show the predicted rest-frame 15 
(j.m LF for t ype-1 AGN and type -2 AGN and compare them with 
the results of lMatute et al.l ( feOOfj) . 

Fig.[8]presents the rest-frame 15 (im LF of type-1 AGN at z = 
0.1 and 1.2. As can be clearly seen, the predi cted results show rea - 
sonable agreements with the measurements of lMatute et al.l ('2006). 
However, it is also clear, especially at z — 1.2, that the predicted 
IRLFs tend to underestimate the number of the most IR-luminous 
type-1 AGN, which is independent of the choices of the evolu- 
tion of HXLF and obscuration. Interestingly, the measurements at 
z — 0.1 can be basically explained by the results predicted from 
most HXL Fs. The only exc eption is the result predicted from the 
HXLF of lAirdetal.1 J20101) modeled with LADE, which signifi- 
cantly underestimated the number of the most IR-Luminous AGN 



even at z — 0.1. Similar to that shown in Fig. [7] these results show 
that the mid-IR LFs predicted from HXLF tend to underestimated 
the number of the most IR-luminous AGN, and become significant 
at 2 > 1. 

Fig. [9] presents the rest-frame 15 um LF of typ e -2 AG N at 
Z — 0.05 and 0.35. As mentioned by IMatute et al] J2006h . the 
observational determination of the rest-frame 15 |i.m LF of type-2 
AGN is much poorer than that of type-1 AGN. So, their measure- 
ment of the density of type-2 AGN can only be considered as a 
lower limit. However, as can be seen in Fig. [9] the results predicted 
from the AGN obscuration evolution mode ls that are constructed 
according to the results of lHasin ger (2008) tend to underestimate 
even the number of type-2 AGN currently measured. Due to the 
much bigger uncertainties in the measurements of the 15 (im LF of 
type-2 AGN, more definitive conclusions cannot be drawn. 



6 DISCUSSION 

By separating the mid-IR LFs of AGN to that for type-1 and type- 
2 AGN respectively, the modeling of AGN SEDs, the evolution of 
LFs and obscuration of AGN can be further constrained. The results 
presented in the Section[5]show that the mid-IR LFs predicted from 
HXLF tend to underestimate the number of the most IR-luminous 
AGN, despite of the general agreements between predictions and 
measurements. This is independent of the choices of the evolution 
models of HXLF and obscuration of AGN, and even more obvious 
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Figure 7. Rest-frame 8.0 |xm LFfor type-1 AGN at z = 1.5, 2.0, and 2.5 as predict ed from the LDDE modeling of HXLF of lUeda et all <2003l) . lEb~rero et alj 
1 2009). lAird et alj J2010 ). and the LADE modeling of HXLF of lAird et alj (2010), respectively. At each panel, the results for six evolution models of fa 
are presented with line styles and colours that are the same as in Fig. [3] The data points show the measured 8.0 \aa LF of type-1 quasars as determined by 
iBrown et al. 1 2006) from a sample consists of 292 24 i^m sources brighter than 1 mJy and selected from Spitzer MIPS survey. The black solid points denote 
the result for objects over the redshift range 1 < z < 5, and the green triangles for those with 1.5 < z < 2.5. All the results predicted from HXLF, especially 
those recently presented, tend to underestimate the number of IR-luminous AGN, independent of the choices of the evolution of HXLF and obscuration. 
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Figure 8. Rest-frame 15 urn LF for type-1 AGN at z = 0.1, and 1.2 as predict ed from the LDDE modeling of HXLF of lUeda et all J20O3h . lEbrero et al J 
J2009ll . lAird et alj l20\(j\, and the LADE modeling of HXLF of lAird et alTj2010l) . respectively. The data points are the measured 15 urn LF of type-1 AGN 
determined bv lMatute et alj J2006I) from a sample of type-1 AGN with redshift in z = [0, 0.2] (top) and z = [0.2, 2.2] (bottom) selected at 15 yon (ISO) and 
12 |j.m (IRAS), and classified by their optical spectra. At each panel, the results for six evolution models of fa are presented with line styles and colours as 
in Fig. [5] The results predicted from all HXLFs t end to underestim ate the number of the most IR-luminous AGN at 2 = 1.2. However, this is not the case at 
Z = 0.1, unless the LADE modeling of HXLF of lAird et alj feOld) is used. 
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Figure 9. Similar to Fig. [8] but for type-2 AGN at z = 0.05, and 0.35. The data points are the measured 15 |xm LF of type-2 AGN determined bv lMatute et alj 
12006) from a sample of type-2 AGN with redshift in z = [0,0.1] (top) and z = [0.1,0.6] (bottom) selected at 15 urn (ISO) and 12 \im (IRAS), and classified 
by their optical spectra. The measurement of the mid-IR LF of type-2 AGN is much poorer than that of type- 1 AGN, and therefore , is mu ch harder to be 
explained. However, the results predicted from the AGN obscuration evolution model constructed according to the results of Hasinger 1 2008) are below even 
the current measurements. 



for the HXLFs recently proposed. Meanwhile, this trend seems not 
significant for AGN with z < 1 and/or less luminous at IR. 

Here, we discuss some possible explanations for this con- 
tradiction between HXLFs and mid-IR LFs. Firstly, this may be 
caused by the missing fraction of AGN, especially those heav- 
ily obscured Compton- thick AGN that cannot be detected by cur- 
rent X-ray observations. Recently. iFu et alj ( 12010m compared their 
mid-infrared spectroscopic selection with other AGN identification 
methods and concluded that only half of the mid-infrared spectro- 
scopically selected AGN were detected in X-ray. However, after 
considering this we find that it only result in a slight improvement 
to the prediction of IRLFs from HXLF Furthermore, this explana- 
tion needs a larger fraction of missing AGN at the high-luminosity 
end, which is in contradiction with the general expectation that 
AGN dominate in the most IR-luminous sources. 

Secondly, the contribution of star formation in AGN host to 
mid-IR, which has not been considered yet, may be important. If 
this is important, the X-ray to mid-IR relation used to predict mid- 
IR LFs from HXLFs in Section [5] needs to be corrected signifi- 
cantly. This is particularly important for sources that are not spa- 
tially r esolved, or with intensive star for mation near the nuclear re- 
gion JLutz et al.l2004l : lHorst et al.l2008T) . We find that if the contri- 
bution of star formation in the host to 8.0 /im and 15 fim emission 
are comparable to the reprocessed nuclear emission, the 8.0 /j,m 
and 15 /im LF predicted from HXLF can be consistent with cor- 
responding mid-IR measurements. Currently, it is still difficult to 
separate the contribution of star formation and AGN to the IR emis- 
sion of galaxies, especially in systems where the two are compara - 
ble and their additive effects are non-linear ( Hopkins et alj|2010h . 
Especially, the relative fractions of their contributions to mid-IR are 
likely to be different in different kinds of galaxies, and may change 
with both luminosity and redshift of the source. However, even for 
powerfully star-formi ng quasars, the co ntribution of star formation 
to mid-IR is small JNetzer et alj|2007l) . So, the possible contribu- 
tion of star formation to mid-IR is not likely the main reason for 
the contradiction. 



Thirdly, the contradiction found in Section [5] may represent 
limitations in the torus model used so far. Although the sim- 
ple torus model inherited from Ballantvne et al. (2006b) has been 
well tested, the distribution and composition of the obscuring 
medium around AGN are still very uncertain. Meanwhile, this 
CLOUDY based torus model essentially assume a smooth dis- 
tribution of dusty obscuring medium. Recently, a clumpy distri- 
bution of dusty obscuring medium is s uggested by some authors 
JNenkova et al.ll2002c iHonig et alJuOOq) . These auth ors have re- 
cently pr oposed sophisticated clumpy torus models qNenkova et alj 
l2008allbl; IHonig et alJ|2O10l ; IHonig & Kishimotoll20ld) that are in 
a better agreement with current IR observations of AGN. Unfor- 
tunately, these clumpy torus models mainly give the IR emission 
properties of AGN, while the self-consistent hard X-ray property is 
not presented. To give the X-ray to mid-IR lumin osity ratios that 
are mo re comparable to the observational results of lMullanev et alj 
( 1201 If) , we have made some imp rovements to the original torus 
model of iBallantvne et alj (2006b). However, the improved torus 
model still have some limitations, which is worth additional efforts 
but is beyond the scope of this paper. 



Finally, as shown in Section [331 the anticorrelation between 
a ox and Ltjv, which has been found by many observations, is im- 
portant for giving X-ray t o mid-IR luminosity ra tios that are more 
comparable to the result of lMullanev et al.U201 lh . Here, we present 
the results obtained by assuming a ox = — 1.5 instead of the typi- 
cal value of —1.4 as used in Section[5] A s can be seen in Fi gs.fTOl 
and 1 1 1 1 the mid-I R LFs measurements of iBrown et al . (2006) and 
iMatute et alj ( 120060 can now be explained much better by the re- 
sults predicted from the HXLFs recently proposed, especially that 
of Ebrero et al. (2009). On the other hand, it is interesting to no- 
tice the dramatical difference betwe en the results in Fig s. [Tol and 
ITTI The 8.0 /im LF measurement of IBrown et al . (2006) is for lu- 
minou s quasars with z > 1, while the 15 fira LF of IMatute et alj 



(2006) is mainly for Seyferts at much lower luminosities and red- 
shifts. While the results at 8.0 fim favor the obscuration evo lution 
models constructed according to the results of lHasinger (2008), the 
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Figure 10. Similar to Fig. [7] but now a ox is assumed to have a smaller value of —1.5 inste ad of the typical val ue of —1.4 as used in Section[5] As can be 
seen, the results have been largely improved after this small change. Now, the measurements of Brown et al. (2006) can be well explained by recently proposed 
HXLF, especially that of Ebrero et al. (2009), when combined with the obscuration evolution models constructed according to the results of Hasinger (2008). 
An even smaller cr ox seems required at vL v (8.0 /im) > 5 * 1O 12 L0. 
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Figure 11. Similar to Fig. [8] but now a ox is assumed to have a smaller value of —1.5 at z = 1.2 instead of the typical value of —1.4 as used in Section 
B] At 2 = 1.2, it is cle ar that the results at highest luminosities have been largely improved. Now, the AGN obscuration evolution models pr oposed by 
Ballan tvne et al] (2006a) seems more favored. However, the change of a ox is not required at z = 0.1, unless the LADE modeling of HXLF of lAird et alj 
l2010l) is used. 



results at 1 5 jim nevertheless give more supports to the models 
proposed bv lBallantvne et alj ( I2006af) . 

These results imply that the obscuration of quasars are dif- 
ferent from that of Seyfert s. Luminous quasars often associate 
with galaxy major mergers dCanalizo & Stockton! 120011) or inter- 



actions JHutchingslll98l iDisnev etalJl 19951 : iBahcall et alii 19971: 
iKirhakos et al.lll9990 . while there are little observational evidences 
for less luminous Seyfert g alaxies being a ssociated with merg ers 
JLaurikainen & Salolll995l ; ISchmitdl200ll ; iGrogin et alj|2005h . If 

the evolution and fueling mechanisms of quasars are very differ- 



© 0000 RAS, MNRAS 000, 000-000 



Evolution ofAGN: connecting X-ray and IR 13 



ent from that of lower luminosity Seyfert galaxies, it is natural to 
expect that the distribution and evolution of the o bscuring medium 
around them are very different. As pointed out bv lBallantvne et al.l 
J2006bh . the dusty mediums obscuring luminous quasars are likely 
distributed in a larger scale and linked to the starburst region, while 
lower luminous quasars and Seyferts are obscured by commonly 
suggested compact torus located at much smaller scale. 



7 SUMMARY 

We have presented a detailed comparison between the 2 — 10 keV 
HXLFs and mid-IR LFs of AGN. The combination of hard X-ray 
and mid-IR provide complementary views for understanding the 
evolution of LFs and obscuration of AGN and their co-evolution 
with galaxies. Four measurements of the HXLF of AGN have been 
collected from the literatures for comparison. A simple but well 
tested torus model, which is based on photoionization and radiative 
transfer code CLOUDY, is then employed to model the composite 
X-ray to IR SEDs for AGN with different luminosities and red- 
shifts. In the modeling of SEDs, we have assumed six evolution 
models of AGN o bscuration, which are constrained by the CXRB 
(Ballantvne et al. 2006a), or constructed according to recent direct 
measurement dHasingen2008l) . The model SEDs of AGN have been 
tested against the observational relations between X-ray and mid - 
IR luminosity of AGN recently given by iMullanev et al.l ( 1201 II) . 
The mid-IR LFs predicted from different combinations of the evo- 
lution models of HXLF and obscuration of AGN are compared 
with the measurements of A GN mid-IR LFs give n bv lBrown et al.l 
J2006h . lMatute et alj J2006I) . and lFuet aljj2010h . respectively. By 
predicting mid-IR LFs for type-1 AGN, type-2 AGN, and total 
AGN from HXLF, and comparing them with corresponding obser- 
vational results respectively, the evolution of LFs and obscuration 
of AGN can be further understood. 

We find that the mid-IR LFs predicted from HXLFs tend to 
underestimate the number of the most IR-luminous AGN, which is 
independent of the evolution model of AGN obscuration. We dis- 
cussed possible explanations for this contradiction. It may partly 
due to the missing fraction of Compton-thick AGN that have been 
predicted by the synthesis model of CXRB, but systematically 
missed by current X-ray observations. However, we find that the 
underestimation to the number of the most IR-luminous AGN can- 
not be eliminated even an extreme assumption, which claims that 
only half of the mid-infrared spectroscopically selected AGN are 
detected in current X-ray observations, is employed. 

We conclude that the contradiction mainly result from lim- 
itations in the modeling of the composite X-ray to IR SEDs of 
AGN. A possible reason is the contribution of star formation in the 
AGN host to mid-IR, which has not been considered yet. We find 
that the contribution of star formation to the 8.0 /jm and 15 /jm 
emission need to be comparable with that of reprocessed nuclear 
emission and even more in the most IR-Luminous sources to elim- 
inate the contradiction. However, the contribution of star forma- 
tion in AGN host to mid-IR is not likely so large but actually de- 
creases with increasing Lir. On the other hand, the contradiction 
may represent limitations in the torus model employed. It is clear 
that the torus model are further constrained to give the specific pre- 
diction of mid-IR LFs for type-1 and type-2 AGN, respectively. 
We have made some im provements to the original torus model of 
iBallantvne et alj d2006bh . such as a different handling of the dif- 
fuse emission and Luv-dependent R ln , to give the X-ray to mid- 
IR luminosity ratios that are more comparable to the observational 



results of IMullanev et al.l ( 120111) . Meanwhile, with some tests we 
find that the anticorrelation between q ox and Luv is important for 
m aking the X-ray to mid -IR luminosity ratios closer to the results 
of Mu llanev et al.ld201lh . Interestingly, a smaller a ox improves the 
prediction of the high-Lm, end of the IRLFs significantly at the 
same time. 

Finally, with all the improvements mentioned above, we find 
that the HXLFs and IRLFs of AGN can be more consistent with 
each other if the obscuration mechanisms of quasars and Seyferts 
are assumed to be different. This is consistent with the idea that 
the obscuration mechanism of luminous quasars dominating at high 
redshifts are very different from that of less luminous Seyferts dom- 
inating at lower redshifts, corresponding to their different triggering 
and fueling mechanisms. However, current measurements of the 
IRLFs of AGN are not accurate enough to allow a more complete 
understanding to this by employing the method presented here. Due 
to this limitation, the conclusions drawn here need to be tested fur- 
ther when better measurements of IRLFs are available. 

More accurate measurements of the IRLFs of AGN, especially 
those determined at smaller redshift bins and more accurately sep- 
arated to that for type-1 and type-2, are very helpful for a more 
complete understanding of the evolution of LFs and obscuration 
of AGN. Based on the observations of newly launched IR space 
telescope such as Spitzer, Herschel and forthcoming James Webb 
Space Telescope (JWST), better measurements of the IRLFs of 
AGN are expected. These measurements will largely improve our 
understanding of the evolution of LFs and obscuration of AGN and 
their co-evolution with galaxies. 
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